Mon. Not. R. Astron. Soc. 000, [[>?? (2011) Printed 2 March 2013 (MN WFjgL style file v2.2) 



Theoretical study of the line profiles of the hydrogen 
perturbed by collisions with protons 



M. G. Santos 1 * and S. O. Kepler 1 

1 Instituto de Fisica, Universidade Federal do Rio Grande do Sul, 91501-900 Porto Alegre, RS, Brazil 



Accepted 2012 January 24. Received 2012 January 17; in original form 2011 September 22 



ABSTRACT 



o 

Ph 

< 

(N 

5? 
00 

Ph 

6 

C3 1 1 INTRODUCTION 



We present theoretical calculations of the quasi-molecular line profiles for the Ly- 
man (Lya, Ly/3, Ly7, LyS) and Balmer (Ha, H/3, H7, H£, He, H8, H9, and H10) series, 
perturbed by collisions with protons. In all calculations we include the dependence of 
the dipole moments as a function of internuclear distance during the collision. The 
broadening from ion collisions must be added to the normal electron Stark broadening. 
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The satellite lines present in the far wing profiles of the Ly- 
man and Balmer series lines of atomic hydr ogen have been 
obser ved in stars and in laboratory plasmas (|Kielkopf et al.l 
2002). Experimentally, the primary difficulties in studying 
these line profiles are (1) the development of techniques to 
observe neutral atomic hydrogen interactions at densities 
high enough for spectral line broadening effects to be signif- 
icant, and (2) quantitatively characterizing the conditions 
under which the spectra are created. On the theoretical side, 
the difficulty arives in developing methods for computing the 
contributions from all possible molecular states of H2 and 
H;J~ to the spectrum of free atoms in collision. 

Satellite features at 1600 A and 1405 A in the 
Lyman— a wing associated with free-free quasi-molecular 
transitions of H2 and are observed in the ultraviolet 
(UV) spectra of a few stars obtained with the International 
Ultraviolet Explorer (IUE) and the Hubble Space Telescope 
(HST). The satellites in the red wing of Lyman—/? are in 
the 905 to 1187 A spectral region, covered by Far Ul- 
traviolet Espectroscopic Explorer (FUSE). The stars which 
show Lyman satellites are DA white dwarfs (WDs), old Hor- 
izontal Branch star s of spectral type A, and A Bootis stars 
jAUard et al.ll2000l ). WDs are the end produc t of evolution 
of ~ 98% of all stars (e.g.. iKepler et al.|[2007l ). Having con- 
sumed all available nuclear fuel (which depends on the indi- 
vidual stellar mass), their source of luminosity is dominated 
by only the residual internal thermal energy. This radiates 
away, exponentially cooling the star, but due to the small 
surface area (a typical WD radius is ~10 9 cm), cooling times 
are long, taking approximately 10 10 years for their effective 
temperatures to decrease from around 10 5 K to near 10 3 K. 
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Therefore, WDs are old, and hidden in their structure is a 
historical record of stellar formation and evolution in our 
Galaxy. 

UV spectroscopic observations of WDs reveal a line 
shape very different from the expected simple Stark broad- 
ening, with satelli te lines in the red w i ng of Lyman— /3 near 
1078 and 1060 A (jKoester et alJll99dll998h . The strengths 
of these satellite features, and indeed the entire shape of 
the wings in the Lyman series, are very sensitive to the de- 
gree of ionization in the stellar atmosphere or laboratory 
plasma, because it is the ionization that determines the rel- 
ative importance of broadening by ion and neutral collisions. 
Therefore the satellites may be used as a temperature diag- 
nostic in a plasma where Saha ionization-equilibrium holds. 
In the case of Lyman— a and the H - H + Lyman—/? satel- 
lites, the shape of the potential plays a dominant role in the 
large difference in th e broadening of quasi-molecular features 
l|Allard et al.lll998bl ). 

Quasi-molecular transitions occur when a photon is ab- 
sorbed or emitted by a hydrogen atom while this atom inter- 
acts with one or more ne ighboring particles, whi ch may be 
atoms, ions or molecules (jAllard fc Kielkopi 1982). At short 
distances, the interaction between the perturbed atoms and 
the other particle(s) usually leads to the forma tion of satel- 
lite mo le cular lines on the wings of atomic lines (jAllard et al.l 
Il998al lbl [l999l ). Structures observed in the wing of Lyman— a 
at 1623 A and 1405 A are due to absorptions of quasi- 
molecular H2 and H^", respectively. The intensity of these 
two satellites depends strongly on the degree of ionization 
in the stellar atmosphere and consequently on the effective 
temp erature T e ff and surface gravity loggr (jHolweger et al.l 
11994) . 

The atmosphere of a WD comprises less than 10 -13 of 
the stellar mass and is dominated by hydrogen (DA white 
dwarfs) or helium (DB/DO) and sometimes shows small 
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amounts of heavy elements (e.g. , DAZ, DBZ and DOZ white 
dwarfs) (see, e.g.. lKoesterll2010f ). The primary technique to 
determine TW and log g of WDs is by comparing observed 
spectra with synthetic spectra from theoretical atmosphere 
models (e.g.. iTremblav et al J [2009V The advantage of this 
spectroscopic technique is that the theoretical line profiles 
are extre mely sensitive to varia tions in the atmospheric pa- 
rameters. iBergeron etafl (|l992h were one of the first to ap- 
ply this to a large number of WDs; they fit the hydrogen 
B aimer lines. The technique was later applied to the UV 
Lyman lines (e.g.. IBergeron et al.lll995l ; iKepler et al.|[l993l ; 
[fl (2005)IV 



Bergeron et"afl (|2002| ) made improvements to the mod- 



els by taking into account the optical quasi-molecular opac- 
ities due to the dipole moments induced by the interaction 
between particles. This interaction causes a temporary re- 
arrangement of electronic charges, which forms a temporary 
molecule, thus producing absorptions or emissions beyond 
those of the isolated atoms. 

Currently there are a large number of optical spectra of 
WDs that have been taken as part of the Sloan Digital Sky 
Surv ey (SDSS) (e.g- lEisenstein et al.ll2006l ; iKleinman et all 
|2010| ). allowing for t he st u dy of the WD mass distri- 
bution (IKepler et all I2007I V IKleinman et all (|2004T l and 
ILiebert et all (120051) derive masse s from sp e ctrosc opic fits, 
and IBergeron et alj l|l99ll) and iKoesterl (|l99ll ) discuss 
the apparent i ncrease in th e aver age mass of DAs with 
T cff <12 000 K. IKepler et alj (120061 ) compare SDSS spectra 
with high S/N spectra obtained with the 8 m Gemini tele- 
scope for four WDs with T e ff ~ 12 000 K and determine that 
the masses derived from fitting the relatively low S/N SDSS 
spectra are systematically overestimated by AM ~ O.13M0. 
Due to a correlation between T e g and log g (a small increase 
in T e ff can be offset by a small decrease in logp), such dis- 
crepancies are concentrated only in the region around the 
maximum of the Balmer lines, 14 000 K > T cB > 12 000 K. 
To stud y the trend of the ap parent increased average mass 
of DAs, IKepler et all (l201Gr ) analyse 1505 such WDs with 
S/N > 20 and T eff > 12 000 K and determined an av- 
erage DA mass of M/M Q = 0.604 ± 0.003. They observe 
that the distribution of the sample is similar to tha t of th e 
Palomar Green sur vey published by ILiebert et al. I (|2005h . 
iFalcon et all |2010l) used an ensemble average of gravita- 
tional redshifts of 449 DAs, which is independent of line 
profil es, and estim a ted a mean mas s of 0.6 4 ± Q.O14M0. 

Koester etafl (120091). IKoesterl (|2010l V ITremblav et al. 
(2009!) and lGianninas et all (|201oF study the possible causes 



for this apparen t increase in mass, not seen in photometric 
deter minations l|Engelbrecht fc Koesterl 120071 ; IKepler et al.l 
|2007| V They conclu de this phenomenon is not caused by He 
contamination, but ITremblav et al] (|201ll ) propose that the 
use of one-dimensional calculations of convection in the at- 
mosphere, compared to more realistic three-dimensional cal- 
culations, might be the cause. There is a huge effort to elimi- 
nate the inconsistencies in the model calculations, and these 
attempts show the need for continued study of the interac- 
tion of hydrogen perturbed by protons, electrons, atoms and 
molecules at high pressure. 

The goal of our work is the study of the interaction of 
the hydrogen atom with protons, which temporarily creates 
the molecule Hj". We calculate the potential of up to 
n = 10, and analyze all possible transitions that can occur, 



as well as the dipole moments for all possible transitions. In 
the literature, we find calculations for only Lya, Ly/3 and 
Ly7 and Ha; for these, we calculate to increased resolution. 
This new theoretical data must be inserted in the modeling 
of line profiles at high density, in addition to the electron 
Stark broadening. 



2 THEORETICAL ANALYSIS 

The dipole moments induced by the interaction between 
radiator and perturber cause a temporary repositioning of 
electronic charges, which temporarily forms a molecule. This 
leads to the possible appearance of satellite features in the 
wing of the line profile and produces different a bsorptions or 
emiss ions than those of non-interacting atoms (|Frommholdl 
[1993(1 . 

The position of the satellite line depends directly on 
the extreme of the potential difference of the transition 
AVext(R) and the dependence of the potential difference 
AV(R) — Vj(R) — Vi(R) with internuclear separation R of 
the states i and j of the allowed transitions of the Hj; the in- 
tensity and shape of the satellites depend on bo th the poten- 
tial d ifference and the radiative dipole moment (|Allard et al.l 
120041 ). 

To calculate spectral line profiles, affected by collisions 
with protons, we must know precisely the theoretical poten- 
tial that describes the molecular interaction between radia- 
tor and perturber, as well as the changes in dipole momenta 
due to the radioactive atom-ion separation for each molec- 
ular state; 

iMadsen fc Peek! (|l97ll ) calculate t he values of the po - 
tential of Hj only for states with n ^ 3- lAHard et afl (|2004l ) 
calculate the potential for states with n ^ 4, which allows 
them to calculate the line profile for Lyman— 7. 

The values of the dipole moments of Hj~ that are avail- 
able in the litera ture were cal c ulated by iRamaker fc Peefl 
l|l972l . ll973l ). and lAllard et al] (|2009l ) refers to the calcula- 
tion of n = 4. 

For this paper, we calculate all potential and dipole mo- 
ments up to n ^ 10 of H^ as a function of the internuclear 
distance R up to 300 a.u. (atomic units), with a step of 0.01 
a.u. We compare the potential of the n ^ 3 states with those 
in the literature (|Madsen fc Peeklll97ll ) and show differences 
smaller than 10 -13 , the precision quoted in the published 
values. We also co mpare the computed dipo le moments with 
those tabulated bvlRamaker fc Peek! (119731 ) and those plot- 
ted in lAllard et al] (2009); they are fully consistent. 



2.1 Potentials for the Hydrogen Molecule Hj 

The structu r e of the molecula r ion Hj~ has been studied by 
iTelleri dl930h.lHvlleraaj Jl93lft . |jafl3 (fl934h.lBa.ber fc Hass* 



(|1935I ). iBates fc Pootsl (|1953| ). and IMadsen fc Peek! (|1971 

but accurate tables of the electronic energy as a function of 
internuclear distance for the excited states n ^ 3 have not 

been published. 

As originally shown by iBurra u (1927), the Schrodingcr 
equation for Hj can be separated and the electronic energy 
E calculated exactly. The Schrodinger equation of the elec- 



tron in atomic units is: 

where r a and ri, are the distances measured in units of ao 
of the electron to the two positive charges a and b, E is 
the energy measured in units of Hartree Eh, and R is the 
internuclear separation. 

In terms of confocal elliptical coordinates, defined by: 
A = (r„ + rb)/R, fi = (r a —rt)/R and 4>, the azimuthal 
angle with respect to the internuclear axis, the Schrodinger 
equation can then be separated in three equations: 
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where *(A, fi, 4>) = A(A)M( M )$(0), p 2 = =^-, e = E-j i 
the electronic energy, and where A and m are the constants 
of separation. 

The solution of equation © is direct, i. e.: 

$ = e (±lm< * ) . (5) 

The solution of equation flgQj, the angular component, is ob- 
tained using the [EM] craa s (1931) method: 



(6) 



where P™+ s (/u) are the generalized Legendre polynomials of 
first order. The summation extends only over even or odd 
values of s, depending on the state calculated. 

Applying equation © in equation ©, we obtain a re- 
currence relation for the coefficients c 3 , 



where: 



b s = 



a s c s -2 + b s c s + d s c s+ 2 = 

2 (s-l)a 

P (2s+2m-3)(2s+2m-l) 

-A - (s + m)(s + m+l)+ 
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The recurrence relation forms an infinite set of ho- 
mogeneous linear equations whose solution requires the de- 
terminant of the coefficient matrix to be null. 

The solution of equation J4]) for the radial part is ob- 
tained using the method of Ijaffel (|l934T l. 



A(A) = (A 2 - l)^(l + A)' T e- pA ^ 
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where a = — — m — 1. 

v 

Applying equation Q to the equation @, the following 
recurrence equation is found: 



where: 

a t = (t+l)(t + m+l), 

/3 t = 2t 2 + (4p - 2a)t - A + p 2 - 2pa - (ra + l)(m + a), 
7t = (t — 1 — <r)(t — 1 — <r — rn). 

The recurrence relation ((9| also produces an infinite set 
of homogeneous linear equations, whose determinant of the 
matrix of coefficients must also be null. When calculating 
the determinant of the infinite recurrence relations of equa- 
tions (0 and ([9]), we obtain two functions fi(A,e,R) and 
f2{A, e, R) which depend on A, e, and R. The electronic en- 
ergy (potential energy Vi) e is obtained by calculating the 
solution of both functions fi(A,e,R) and f%(A, e,R), e.g., 
by using the Newton-Raphson method. Practical consider- 
ations require the truncation of the recurrence relations (JSJ 
and ©. We choose t=s=25 for the data presented here, be- 
cause larger values for t and s do not produce any significant 
change in A and e. The coefficients c s and gt are calculated 
using the method of continued fractions. 



2.2 Electric dipole moments 

Choosing the origin of the coordinate system in the center 
of the line segment which unites the two positive charges 
and taking this line as the axis z, the matrix element to be 
evaluated of the components of the electric dipole moment, 
Dij, specified by: 



Dij =e(i\r\j) = eQij{r) , 



(10) 



where e represents the electron charge, r is the position vec- 
tor of the electron relative to molecular fixed frame, in co- 
ordinates x, y and z of the electron, and tyj are the elec- 
tronic functions in the Born-Oppenheimer approximation of 
the states i and j and Q;,j(r) the matrix elements. 

The matrix elements Qij (r) of the va rious components 
of the electric dipole moment according to lRamaker fc Peek! 
l|l973h . are specified by: 

Qi,j (r) = J ** (A, [i, 0) (r) tf, (A, p, 0)dr . (11) 

where r is the volume element. Qi,j{r) is obtained by insert- 
ing the equations ((5j), ((6]l and (JHJ into equation (|lip . The 
matrix elements Qij(r) can be written as a function of the 
coordinates x, y and z as follows: 



Qi,j(x) = NiNj 
QiAv) = N t Nj 
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atgt+i - Ptgt + itgt-i = 0. 
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where ilij (a, b, c, d) is defined by: 



n <J -(a I &,c,d) = / M*(m)(1 - fJ, 2 ) a (J. b M j (n)dn 
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The normalization constants can be calculated from: 



Ni = ( y ) [a,<(0, 0, 0, 2) - Q M (0, 2, 0, 0)] 
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After calculating the values of Q,,j(r ) for each compo- 
nents x, y and z, we obtain .Djj as follows l|Herman fc Wallia 
Il956h : 



\Di 



\eQij(x)\ 2 + \eQ^(y)\ 2 + \eQ l , J (z)\ 2 



(18) 



We compute the values of the electric dipole moments, 
equation (|18[1 . for all transitions of the Lyman and Balmer 
series for the H^" molecular ion for internuclear separations 
up to 300 a.u. . 

Table [1] shows values of the el ectric dipole moments o f 
this work compared with those of iRamaker fc Peek! ()l973l ) 
for some transitions, demonstrating that our work presents 
a higher precision and greater coverage of the internuclear 
distance. 

For equations (|6]) and (£Sj) to be solutions of the ap- 
propriate electronic functions tyi, the summations must be 
allowed to go to infinity. However, we achieve convergence 
of the series for s = 31 and t = 10. 

The data of the potential and electric dipole moments 
are plotted in subsequent sections for the Lyman series (Lya 
to Ly5) and Balmer-a. Electric dipole moments are plotted 
for Balmer-a and Balmer-/?. 



3 PROFILES AND SATELLITES 

Currently there are different approaches to the calculation of 
the line profiles, namely, quasi-static, quantum-mechanical 
and unified theory. 

In the quasi-static approxima tion, the expression of the 
cross-section may be written as (iMargenau fc Lewis! 1 19591 ; 
lAllard fc KielkopJll982l ; I Rohrmann et al.ll201ll ) 



16 3 

— tv n v a 

3 P 



hv 



\dVij/dR\ 



(19) 

where ft — p=, k is the Boltzmann constant, and T is the 
gas temperature. n p is the mean density of perturbers in 
the gas, Vij the potential energy difference associated with 
a transition with an energy hv = Vij = Vj(R) — Vi(R) at 
internuclear separation R, Dij(R) is the dipole moment, and 
Ei(R) corresponds to the change in energy (with respect to 
R — oo) in the state i of the Hj. When there is an extreme 
in the curve corresponding to the difference in the potential 



Table 1. Tab l e of comparison between the data of 
IRamaker fc Peekl J 197^ t and this work! of the electric dipole 
moments for Lyman-Q of the HjJ" in atomic units. 



_n,^UQ ) 


l V !^ 11 O 1 \ 1 C\1~\ 

_L 1 £LI1£>1 LHJ11 






1.0 


2pn u -lsag 


0.5494 


0.5493774 


11.0 


2pir u -lsa g 


0.7169 


0.7168555 


60.0 


2pir u -lsa g 


0.7450 


0.7449931 


1.0 


3pa u -lsa g 


0.1651 


0.1651107 


11.0 


3po u -lsa g 


0.4859 


0.4858565 


80.0 


3pcr u -lscrg 


0.5207 


0.5207119 


1.0 


4f<r u -lsa g 


0.001055 


0.0010541 


11.0 


Afa u -lsa g 


0.6609 


0.6608712 


80.0 


4/fTu-lSCTg 


0.5326 


0.5325508 


1.0 


2so g -2po u 


1.241 


1.2406089 


11.0 


2sa g -2p<T u 


0.4844 


0.4844248 


80.0 


2sa g -2pa u 


0.5207 


0.5207119 


1.0 


3dcj g -2pa u 


1.010 


1.0104411 


11.0 


3d<jg-2pa u 


0.6349 


0.6349347 


80.0 


3da g -2pa u 


0.5326 


0.5325508 


1.0 


3dn g -2pa u 


0.8584 


0.8584398 


11.0 


3dn g -2pa u 


0.7935 


0.7935359 


60.0 


3diTg-2pcr u 


0.7450 


0.7449931 



energies between the upper and lower states, is zero. 

Therefore, if the derivative is zero, there is a divergence, 
called a satellite line. 

In the quantum-mechanical approach, the cross section 
for the s pontaneous plus stimulated association process is 
given by (|Zygelman fc Dalgarnoiri990l ) 



64 7T J 



3 c 3 k? 1 - exp{-hv/kT) 



xp[N M 



+ (N +1)M 2 ,, 



(20) 



where E is the relative collision energy, k is the wave num- 
ber of relative motion, p is the probability of approach in 
the initial electronic state, N is the initial rotational quan- 
tum number, v is the final vibrational quantum number, 
and M is the electric dipole matrix element connecting the 
initial continuum and final rotational-vibrational states of 

the ground electronic state. 

The unified theory has been developed in lAllard et al.l 

(1999), and a detailed discussion is presented there. The fun- 
damental expression of the normalized sp ectrum line F v (Av ) 
is computed from the Fourier transform (|Allard et alj|2009h 



F„(Av) = -Re 

7T 



$(s)e 



s ds. 



(21) 



Figure [T] shows the comparison of the line profile of 
Lya of the three approaches mentioned above. The solid 
line shows the approach from the unified theory, which has 
a resolution determined by Fourier transform (equation [21]), 
while the quantum-mechanical approaches (dashed line) and 
quasi-static (solid line) values provide similar but distinct 
profiles due to the calculation of the derivative by the numer- 
ical difference in the quasi-static and by the series expansion 
in quantum-mechanics. The three profiles are similar; how- 
ever, the quasi-static approach better represents the correct 




global shape of the line wings and position of t he satellites, 
in co mparison with previously published data (|Allard et al.l 
1999). As the main aim of our work are: the calculations 
of the electric dipole moments, the potentials of H^" ion, 
the calculations of the global shape of line profiles and the 
location of the satellites, we use the quasi-static approach 
to calculate the line profiles. Unfortunately, the quasi-static 
approach does not provide an appropriate shape of the cen- 
ter of the lines, thus there is a need for better modeling of 
the center lines of each transition. 

Our interest is to find the correct position of the satel- 
lites in the wings of the line profile and their general form 
(i.e., the divergences in the profile) and bring these results 
to the attention of astrophysicists. If they are observable, 
they will be specific for protons (although other ions may 
give qualitatively similar effects). 

The method we u sed to obtain the profiles following 
iRohrmann et al.l l|201lf ). For lines profile of the Lyman and 
Balmer series in a quasi-static approach, we identified each 
pair of levels of the quasi-molecule H - H + that con- 

tributes to the formation of a given line. Only transitions 
between states u (ungerade) and g (gerade) are allowed, i.e., 
with opposite parity and taking m — m! |= or 1 as the 
allowed dipole. The second step is to calculate the effective 
cross section aij{v), equation (|19[) using the potentials and 
dipole moments for the transitions already identified. For 
the calculation of (Jij{v) we need to analyze the frequen- 
cies v that arise from the difference in energy between these 
levels hv — Vj — Vi. However, we must take into account 
that different internuclear distances can produce the same 
frequency. The effective cross section is the sum of each of 
these contributions. 

The line profile is the sum of the effective cross sections 
(Jij(v) multiplied by the statistical weights of each of the 
respective transitions. The calculated profile was convolved 
with a Maxwellian distribution of velocities to take into ac- 
count the thermal Doppler broadening. 

The line profiles and satellite lines of Lyman and Balmer 
series shown in this paper are for single densities of 10 to 
10 18 cm -3 , typical of the atmosphere of white dwarf stars 
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Figure 3. Lyman-a profile with density of H + perturbers (T= 
12 000 K), n+ = 10 16 , n+ = 10 17 and n+ = 10 18 cm" 3 (from 
bottom to top). 

and that can be produced in z-pinch laboratories. We remind 
the reader that the real profile must include the electrons 
Stark broadening, which do not have satellites. 

3.1 Lyman series 

3.1.1 Lyman-a 

The Lyman-a total profile depends on 6 individual transi- 
tions. The line profile calculations shown in Figures [2] and [3] 
have been done for temperatures of 4 000 K and 12 000 K, 
at proton densities of 10 16 , 10 lr and 10 18 cm" 3 . 

Figure [4] shows the potential energy for the transitions 
that contribute to the profile of Lyman-a and Figure[5]shows 
the electronic transition dipole moments for Lyman-a due 
to H perturbed by collisions with protons. 

Table [2] lists the wavelengths of the satellite lines with 
the upper and lower state identifications and the distance 
of the atom-ion pair at which the potential extreme oc- 
curs. Several of these features can be identified in the total 
Lyman-a profile of Figures [2] and [3] 



Table 2. Allowed transitions and satellites due to H - H + colli- 
sions of Lyman-a to the distance R(A) of the atom-ion pair at 
which the potential extreme occurs. 



Label 


Upper 
level 


Lower 
level 


(cm 1 ) 


(A) 


R 

(A) 


1 


2pw u 


lsa g 


-1628.54 


1240.56 


4.72 


2 


3pcr u 


ISCTg 








3 


4/o- u 


lsa g 


-1172.40 


1233.58 


11.07 


•1 


2sa g 


2pa u 


8139.91 


1106.47 


2.65 


5 


3d<Tg 


2pa u 


-11055.90 


1404.87 


4.56 








16001.33 


1017.93 


1.39 


6 


3dir g 


2pa u 


21384.28 


965.05 


1.67 




Figure 4. Transitions difference potential as a function of inter- 
atomic distance of the H^" of Lyman-o. 



We l ist in Table [2] the sat ellite around A = 1400 A, ob - 
served bv lKoester et all (|l985l ) and lNelan fc Wegnerl (|l985l ). 
who proposed the lines are caused by quasi-molecular ab- 
sorption of the molecule. 
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Figure 5. Electronic transition dipole moments of Lyman-a due 
to H perturbed by collisions with protons. The D(R) is given here 
in atomic units. 
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Figure 6. Lyman-/3 profile for T= 4 000 K and density of H + 
perturbers of = 10 16 , = 10 17 and = 10 18 cm -3 (from 
bottom to top). 
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Figure 7. Lyman-/3 profile for T= 12 000 K, and density of H + 
perturbers of = 10 16 , ritj = 10 17 and = 10 18 cm -3 (from 
bottom to top). 



3.1.2 Lyman-f3 

The profile of Lyman-/3 depends on 10 transitions. Table [3] 
lists all the wavelengths of these transitions and the pos- 
sible lines satellites. Satellites of Lyman-/3 were observed 
in the spectra of the DA white dwar f Wolf 1346 with th e 
Hopkins Ultraviolet Telescope (HUT) (|Koester et al.lll996h . 
iKoester et al.l (|l998h observed satellites absorption features 
at A = 1060 A and A = 1080 Ain the spectra of 4 other 
targets obtained with the Orbiting and Retrievable Far and 
Extreme Ultraviolet Spectrometer (ORFEUS). These satel- 
lites are due to the contributions of the transitions 4fn u - 
lsa g and 5ga g -2pa u of the profile of Lyman-/3, a ccording to 
the profiles of calculated bv lAllard et alj (Il998al ). 

The Ly/3 line profiles shown in Figuress [6] and [7] are 
for temperatures of 4 000 K and 12 000 K, and densities of 
protons 10 16 , 10 17 and 10 18 cm" 3 . Figure [8] shows all the 
differences of the potential energy of the transitions that 
contribute to the Lyman-/? line profile. In Figure [5] we plot 
the dipole D(R) for transitions of the Lyman-/? due to H 
perturbed by collisions with protons. 



Table 3. Allowed transitions and satellites due to H - H + col- 
lisions of Lyman-/3 to the distance R(A) of the atom-ion pair at 
which the potential extreme occurs. 



Table 4. Allowed transitions and satellites due to H - H + col- 
lisions of Lyman-7 to the distance R(A) of the atom-ion pair at 
which the potential extreme occurs. 



Label 


Upper 


Lower 






R 




1 vel 


level 


1 

(cm ) 




1 41 


1 


3pn u 


ISOg 








2 


4pa u 


lsa g 








3 




ISCTg 


-3361.25 


1062.64 


9.85 


4 


5fa u 


lsa g 








5 


Qha u 


lsa g 


-1037.33 


1037.03 


21.45 


6 


3sa g 


2pa u 


22200.04 


835.65 


1.96 


7 


Ada g 


2pa u 


-633.34 


1032.71 


9.46 








25868.59 


810.80 


1.55 


8 


MlTg 


2pa u 


28666.11 


792.82 


1.70 


9 


5ga g 


2pa u 


-4899.9272 


1080.31 


12.65 








38880.59 


733.42 


1.69 


10 


5gir g 


2pa u 


-516.61 


1031.46 


18.89 








38914.86 


733.24 


1.69 
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Figure 8. Transition potential energies as a function of inter- 
atomic distance of the H~j" of Lyman-/3. 
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Figure 9. Electronic transition dipole moments of Lyman-/3 due 
to H perturbed by collisions with protons. The D(R) is given here 
in atomic units. 



Label 


Upper 
level 


Lower 
level 


l 

(cm ) 


H J 


R 

( 4 \ 


1 




ISOg 




- 


- 


2 


bpa u 


lsa g 




_ 


_ 


3 


5/7r u 


lsa g 


- 


974.62 


16.65 


4 


6fa u 


lsa g 


- 


- 


- 


5 


Qhn u 


lsa g 


-2188.78 


993.96 


20.69 


6 


7ha u 


lsa g 


-230.90 


974.98 


29.66 


7 


8ja u 


lsa g 


-784.17 


980.27 


36.06 


8 


4sa g 


2pa u 


30153.49 


752.16 


1.82 


9 


5da g 


2pa u 


32004.30 


741.83 


1.62 


10 


5dn g 


2pa u 


33520.98 


733.58 


1.70 


11 


6ga g 


2pa u 


-1395.38 


986.18 


17.79 








39005.56 


705.21 


1.69 


12 


Gg-n-g 


2pa u 


39025.78 


705.11 


1.69 


13 


1ia g 


2pa u 


-2627.07 


998.31 


25.06 








42340.74 


689.00 


1.70 


14 


7in g 


2pa u 


-565.52 


978.18 


31.65 








42342.54 


688.99 


1.70 



3.1.3 Lyman-^f 

There are 14 allowed transitions due to collisions of H - H + 
and the satellites predicted to occur are listed in Table |U 
The satellite absorption features at A = 995 A was observed 
in a HUT spectrum and later in spectra recorded by OR- 
FEUS. This Lyman-7 satellite was also visible in the spec - 
tra of the wh ite dwarf CD - 38° 1 0980 ll Wolff et aLllgOOlT) 
and Sirius B l|Holberg et al. l l2003h . lAllard et all (|20o3 7 cal- 
culated the profiles of Lyman-7 une based on accurate ab 
initio potentials to explain the observation of this feature. 

The line profile calculations shown in Figure [10] and 
Figure [11] have been evaluated at a temperature of 4 000 K 
and 12 000 K for differents densities of protons 10 16 , 10 lr 
and 10 18 cm" 3 . 

Figure [T5] shows all the potential energy differences for 
the transitions that contribute to the Lyman-7 une profile 
and Figure[l3]plots the electronic transition dipole moments 
for Lyman-7 due to H perturbed by collisions with protons. 



3.1.4 Lyman-S 

There are 18 allowed transitions due to collisions of H and 
H + , and the satellites predicted to occur are listed in TableJS] 

Figures [141 and [T31 show the line profiles as a function of 
temperature and density of perturbers. 

Figure [TS] shows all the potential energy differences for 
the of the transitions that contribute to the line profile of 
Lyman-5. Figure [17] shows the electronic transition dipole 
moments for Lyman-<5 of the and Figure [TS] shows the 
total line profiles of Lya, Ly/3, Ly7 and LyS perturbed by 
protons for a single T= 12 000 K and n+ = 10 17 cm~ 3 
This f igure is similar to Figure 3 obtained bv lHebrard et all 
|201ll) . showing that our calculations and approximations 
are consistent with data already published. 
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Figure 11. Variation of the Lyman-7 profile with density of H+ 
perturbers (T= 12 000 K), n+ = 10 16 , n+ = 10 17 and n+ = 10 18 
cm -3 (from bottom to top). 




Figure 12. Transition potential energy differences as a function 
of interatomic distance of the Hjj" of Lyman-7. 



Table 5. Allowed transitions and satellites due to H - H + col- 
lisions of Lyman-5 to the distance R(A) of the atom-ion pair at 
which the potential extreme occurs. 



Label 


Upper 


Lower 




v 


R 




level 


level 


[cm 1 ) 


(A) 


(A) 


1 


5p7r u 


XatTg 








2 




XstTg 








3 


6/7T„ 


XstTg 








4 


7/tT„ 


lSCTg 








5 


7hir u 


XSUg 


-679.23 


956.16 


27.56 


6 


8hcr u 


XatTg 








7 


SjTVu 


XSCTg 


-1413.57 


962.93 


36.17 


8 




XStTg 


-312.56 


952.82 


44.76 


9 


X0la u 


XSUg 


-582.87 


955.28 


54.98 


10 


5SCTg 


2pa u 








11 


6da g 


2po u 








12 


ddiTg 


2pa u 








13 


7go g 


2pa u 


-76.95 


950.69 


26.02 


11 


Tg-Kg 


2ptr u 








15 


8i(Tg 


2pa u 


-1119.44 


960.21 


31.54 


16 


8in g 


2po u 


-99.84 


950.90 


42.38 


17 


9k<Tg 


2pcr u 


-1589.40 


964.56 


41.91 


18 


9klTg 


2pa u 


-475.45 


954.31 


49.00 



3.2 Balmer series 

The Balmer lines have transitions that must obey the same 
golden rules as the Lyman series transitions. 

3.2.1 Balmer-ct 

The total profile of Balmer-Q depends on the 32 individual 
transitions. The line profile calculations shown in Figure [23] 
have been calculated at a single temperature of 12 000 K 
and protons density of 10 17 cm" 3 . 

Table [S] lists only the transitions that generate possi- 
ble satellite lines on the wings of the Balmer line profiles. 
In Figure [19] we plot the transition dipole moments for the 
satellites listed in Table [Td] In laboratory observations, mea- 
surements of the spectrum of a laser-produced plasma con- 
firm that the A = 8650 A satellite exists in the Ha red 
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Figure 14. Lyman-<5 profile for density of H + perturbers (T= 
4000 K), n+ = 10 16 , n+ = 10 17 and n+ = 10 18 cm" 3 (from 
bottom to top). 
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Figure 17. Electronic transition dipole moments of Lyman-<5 due 
to H perturbed by collisions with protons. The D(R) is given here 
in atomic units. 
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Figure 15. Lyman-<5 profile for density of H+ perturbers (T= 
12 000 K), n+ = 10 16 , n+ = 10 17 and n+ = 10 18 cm" 3 (from 
bottom to top). 
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Figure 16. Potential energy differences for the allowed transi- 
tions, as a function of interatomic distance of the H J of Lyman-<5. 




Figure 18. Total profile of Lyman-a, Lyman-/?, Lyman-7 and 
Lyman-<5 perturbed by protons (n~^ = 10 17 cm -3 ) at a tempera- 
ture T= 12000 K. 



wins (|Kielkopf et al.ll2002l ) due the transitions Sgcg-AfcTu 
and 4fn u -3diTg. 

Figure 1201 shows the comparison between the VCS 
l|Vidal et alj Il973l ) profile calculated by Koester (private 
communication) for T= 12 000 K and density 10 16 cm -3 , 
with the profile of calculated for the same temperature 
and density. It is observed that at high energies both profiles 
coincide asymptotically, while, due to possible formation of 
low-energy molecules, discrepancies arise between the pro- 
files. 



3.2.2 Balmer-13 

The total profile of Balmer-/? depends on the 46 individual 
transitions. Figure [2T1 shows the profile of H/3 calculated for 
protons density of 9.324 x 10 16 cm" 3 and 1.554 x 10 18 cm" 3 
with T — 13 526 K. We can see the center of the line is split 
and asymmetric. 

Table [7] lists the lines satellites due to the H/3 and Fig- 
ure [221 plots the respective transitions listed in the same 
table. The line profile calculations shown in Figure [23] is for 
temperature of 12 000 K and protons density of 10 17 cm" 3 . 



Table 6. Satellites due to H - H + collisions of Balmer-a to the 
distance R(A) of the atom-ion pair at which the potential extreme 
occurs. 



Label Upper Lower Av + A + R 

level level (cm -1 ) (A) (A) 



1 


5ga g 


4f°u 


-3844.18 


8783.59 


13.10 








4956.58 


4954.01 


5.60 


2 


5gcr g 


2pn u 


-4996.57 


9772.81 


12.66 


3 


6ha u 


3dcr g 


40925.21 


1780.80 


3.17 








-702.54 


6883.97 


23.03 


1 




3d(Tg 


-1271.27 


7164.46 


13.13 








23797.73 


2562.34 


3.06 


5 


Ad<jg 


3p<Ju 


-3214.16 


8323 01 


6.47 








12145.77 


3652.99 


1.01 


6 




3dw g 


-3658.85 


8642.89 


9.33 








7163.74 


4465.72 


1.69 


7 


if-Ku 


2S(Tg 


-5245.35 


10016.33 


8.69 








37537.98 


1895.12 


1.00 


8 


5f(T u 


2S(Tg 


-424.84 


6754.83 


15.47 


9 


Gha u 


2S(Tg 


-1500.75 


7284.22 


20.33 


10 


AdtTg 


2pit u 


-634.82 


6852.02 


10.20 


11 


bg-rtg 


2pir u 


-601.70 


6836.51 


18.74 


12 


5ga g 


3pa u 


-6042.60 


10885.60 


12.08 








23287.20 


2596.30 


1.18 


13 


5gw g 


3pa u 


-1130.91 


7093.14 


16.58 








23303.15 


2595.23 


1.18 


14 


Qha u 


3dir g 


-1117.22 


7086.25 


21.38 








23163.48 


2604.67 


2.10 




Figure 19. Transition dipole moments for the satellites of 
Balmer-a due to H perturbed by collision with protons. The tran- 
sitions plotted here are listed in Table [6] The D(R) is given here 
in atomic units. 



3.2.3 Balrner-'y 

The total profile of Balmer— 7 depends on the 60 individ- 
ual transitions. The line profile shown in Figure [23] is for a 
temperature of 12 000 K and proton density 10 17 cm" 3 . 



3.2.4 Balmer-5 

The total profile of Balmer-5 depends on the 74 individual 
transitions. The line profile calculations of Figure [23] have 
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Figure 20. Comparison between the profile of Balmer-a of HJ 
with the profile VCS - electron Stark broadening, for T= 12 000 K 
and density 10 le cm" 3 . 




9.324 X 10 16 cm -3 of the perturbers and T = 13 526 K. Dashed 
line is the Balmer-/3 profile with density = 1.554 X 10 18 cm -3 
of the perturbers and T = 13 526 K. The line profile can be used 
for diagnostic of the physical conditions of a laboratory measure- 
ments. 



been done at a temperature of 12 000 K for the protons 
density 10 17 cm" 3 . 



3.2.5 Balmer-t 

The total profile of Balmer-e depends on the 88 individual 
transitions. The line profile calculations of Figure [23] have 
been done at a temperature of 12 000 K for the protons 
density 10 17 cm" 3 . 



3.2.6 Balmer-8 

The total profile of Balmer-8 depends on the 102 individual 
transitions. The line profile calculations of Figure [22] have 
been done at a temperature of 12 000 K for the protons 
density 10 17 cm" 3 . 



Table 7. Satellites due to H - H + collisions of Balmer-/? to the 
distance R(A) of the atom-ion pair at which the potential extreme 
occurs. 





Upper 


Lower 






n 




level 


level 


/ — 1 \ 

(cm A ) 


(A) 


(A) 


1 


5/ 7T-u 


1S(Jn 


-992.00 


5110.59 


13.83 


2 


6/l7T- u 


a 


-2674.11 


5591.24 


19.90 


3 


5g3g 


2p7T u 


-1466.46 


5237.59 


16.80 


1 




2p7T u 


-1483.16 


5242.17 


17.74 


5 






-2040.74 


5400.01 


16.50 








23365.48 


2276.62 


1.19 


6 


7i(Jn 

i v\j g 


3p<r u 


-2977.17 


5687.62 


24.59 








26648.74 


2118.28 


1.19 


7 


7 ill n 


3p(T u 


-816.65 


5065.20 


29.99 








26649.62 


2118.24 


1.19 


8 




3d(Jg 


-1811.98 


5334.12 


21.66 








35621.24 


1779.97 


3.17 


9 


8j<T u 


3d(Tg 


-381.85 


4956.05 


29.99 








41130.19 


1621.02 


3.19 


10 


5f7Vu 


3d,7Tg 


-285.48 


4932.49 


16.46 








12088.84 


3062.96 


1.91 


11 






-2269 96 


5467 69 


20.64 








17843.96 


2603.95 


2.11 


12 


6hS u 


3d7Tg 


-308.92 


4938.20 


28.53 








17875.84 


2601.79 


2.13 


13 


7hcr u 


SdlTg 


-302.86 


4936.72 


29.57 








21106.24 


2400.06 


2.12 


11 


Sj(T u 


3dlVg 


-853.70 


5074.72 


36.03 








23255.19 


2282.35 


2.15 


15 


5da g 




-11612.48 


11177.27 


1.95 








3925.46 


4084.18 


8.17 


16 


6ga g 




-878.51 


5081.12 


20.00 








6602.41 


3681.66 


6.08 


17 


7iTT g 




-405.68 


4961.91 


29.99 








12418.51 


3032.34 


6.81 
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Figure 22. Electronic Transition dipole moments for the satel- 
lites of Balmer-/? due to H perturbed by collision with protons. 
The transitions plotted here are listed in Table 171 The D(R) is 
given here in atomic units. 
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Figure 23. Profiles of the Balmer series (Ho, H/3, H7, H5, He, H8, 
H9, and H10) for ion density n+ = 10 17 cm" 3 and T= 12 000 K. 



3.2.7 Balmer-9 

The total profile of Balmer-9 depends on the 116 individual 
transitions. The line profile calculations of Figure [23] have 
been done at a temperature of 12 000 K for the protons 
density 10 17 cm -3 . 



3.2.8 Balmer-W 

The total profile of Balmer- 10 depends on the 130 individual 
transitions. The line profile calculations of Figure [23] have 
been done at a temperature of 12 000 K for the protons 
density 10 17 cm" 3 . 



4 DISCUSSION AND CONCLUSION 

Satisfactory theory and data for the line profiles do exist 
due to the el e ctron Stark broad ening of neutral hydrogen 
jLemkd Il997l ; IVidal et all ll973Ti, and f o r 21 optical lines 
of ne utral helium ( Barnard et al.l 1 19691 ; iBeauchamp et all 
1997). The first three Lyman lines of H broadened by ion- 
ized and neutral perturbers, including a number of satellite 
features, are well described by the work of Nicole Allard 
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Figure 24. Total profile of the Balmer for T= 12 000 K and 
n+ = 10 17 cm" 3 . 



and collaborators (e.g. lAllard et al. I l2009l . and papers cited 
therein), and the atomic data are obtained from a number 
of atomic databases, pr edominantly the line lists from Ku- 
rucz and collaborators (|Kurucz fc Bell 1 19951). and the Vi- 
enna Atomic Line Database (VALD) llKupka et al 1 1 19991 . 
l2000l ; lRvabchikova et al.lll997l ; |Piskunov et al.lll995h . For all 
other cases the situation is much less satisfactory, so there 
is a need for theories and data to explain the line profiles 
and satellite lines. 

We extended the potentials dipole and line profile cal- 
culations to n=10 for the collision of H on H + perturbers 
important at high densities, as observed in the atmosphere 
of white dwarf stars. 

The laboratory observations of a laser-produced plasma 
confirm that satellites appear on Lyman— a due to collisions 
with neutral atoms and protons. This experimental confir- 
mation of the theory supports the identification of these fea- 
tures in the Lyman— a spectra of white dwarf and A Bootis 
stars. In high temperature plasmas such as laser-produced 
plasmas, dielectronic recombination process is very impor- 
tant. This process strongly affects the ion abundances and 
produces dielectronic satellite lines. The satellite lines ap- 
pear near the resonance line through spontaneous radia- 
tive decay from the autoionizing states. Satellite line emis- 
sions from highly ionized ions are often used for plasma 
diagnostics which are important in various kinds of plas- 
mas, e.g., To kamak, inertial fusion, a strophysical and plasma 
sources, etc (|Yamamoto et ai1l2004l ). The comparison of the 
Lyman— a wing with line shape models is a tool for deter- 
mining neutral and proton densities in a hydrogenic plasma. 
These experiments also confirm that the variation of the ra- 
diative dipole moment is an important factor in determining 
the far wing emission of Lyman— a. When D(R) differs sig- 
nificantly from its asymptotic value at an R close to the 
region forming a satellite, the strength of the wing may be 
enhanced (or diminished) considerably. 

In a collision-induced spectrum such as this one, the 
dipole moment may be more important than the potential 
in determining the shape of the satellite. These many-body 
perturbations alter the far line wing, adding multiple satel- 
lites and producing a strong continuum from the vacuum 
ultraviolet to the ultraviolet and visible. 



The present calculations are done in an adiabatic ap- 
proximation using a rectilinear trajectory. This should affect 
slightly the shape of the satellite, although no significant dif- 
ferences are expected. 

The sum of the profiles calculated of Lya, Ly/3, Ly7 and 
Ly<5, are shown in Figure [TH] and the sum of the profiles of 
Balmer series in Figure [24] to T= 12 000 K and n+ = 10 17 
cm" 3 . The ion collision profiles differs significantly from the 
electron Stark broadening, showing not only satellite lines 
but also significant wing values that affect the opacity. 

The data is available in our web page at 
http:/ /astro. if.ufrgs.br/marcios 
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